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Abstract 

This report describes a method for the primary calibration of phase 
noise and amplitude noise measurement systems. 

In the field of metrology, the term primary refers to a standard whose 
quantity value and measurement uncertainty are established without re- 
lation to another measurement standard for a quantity of the same kind; 
or to a procedure used to realize the definition of a measurement unit 
and obtain the quantity value and measurement uncertainty of a primary 
measurement standard. Accordingly, a phase and a normalized amplitude 
have to rely on ratio measurements. 

This is a working version, which I expect to update in the following 
six months adding some crucial experimental results. 
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1 Introduction 

As physics and technology progress, phase noise becomes an increasingly rel- 
evant issue in a number of fields. The demand for lower phase noise is seen 
in a number of examples, like high speed electronics, and telecommunications 
Extreme low phase noise specifications are found in military and space oscilla- 
tors. Some types of military radars can achieve the desired resolution only if the 
phase noise of the X band reference signal does not exceed 1CP 14 rad 2 /Hz (—140 
dB) at 10 kHz off the carrier. Space research requires a frequency stability in 
the lower 10~ 15 at r = 1 . . . 10 3 s, which means a phase noise exceed 10~ 13 
rad 2 /Hz (—130 dB) at 1 Hz off the carrier if the signal is taken at the standard 
100 MHz output. Of course, we expect that the need for measurement accuracy 
will inevitably follow the specifications for low phase noise, and consequently for 
high measurement sensitivity. Besides, primary metrology is around the corner. 
In fact, in 2010 the definition of the second is expected to switch 1 from the 
F = 3 <-> F = 4, Am/ = hyperfine transition of the 133 Cs atom (9.192 631770 
GHz) to an optical transition, to be determined. In this business, we expect 
that the phase noise of some synthesizers and transfer oscillators is a part of the 
metrological chain, for it has to be understood with high accuracy. 

This article privileges phase noise vs. amplitude noise because of the needs 
and of the culture of our laboratory. Nonetheless, the method proposed here 
relies upon amplitude noise measurements, and provides a solution in which AM 
and PM calibration have the same accuracy. 

In the field of metrology, the term primary [28, 8, 15] refers to a standard 
whose quantity value and measurement uncertainty are established without re- 
lation to another measurement standard for a quantity of the same kind; or to 
a procedure used to realize the definition of a measurement unit and obtain the 
quantity value and measurement uncertainty of a primary measurement stan- 
dard. By contrast, secondary means that the standard, or to the measurement 
procedure, relies upon calibration against, or comparison with, a primary mea- 
surement standard for a quantity of the same kind. It is worth mentioning that 
an empirical unit can have higher reproducibility than the primary counter- 
part, as it happens with the Josephson and with the quantum-Hall electrical 
standards, yet inferior absolute accuracy as related to the SI. 

1.1 Basic definitions 

The SI unit of angle, the radian (rad) is now considered a derived unit because an 
angle can always be defined in terms of the ratio of two homogeneous quantities. 
Formerly, it was considered an auxiliary unit. In the domain of electric circuits, 
it is usual to represent the sinusoidal signal v(t) = \/2V rms cos[u>oi + <p] as the 
complex number V — V^^e^ obtained by freezing the frequency luq. This 
complex number is called 'Fresnel vector', 'phase vector', or 'phasor'. When V lms 
and ip are allowed to vary (slowly), the term 'analytic signal' is used instead. 
The object of our interest is the primary measurement of the phase <p, focusing 
on small randomly fluctuating angles. 

lr This was said at the panel discussion on the future of SI at the CPEM conference, Torino, 
Italy, July 9-14 2006. Yet, it seems that this decision is still not in the official documents 
(thanks to Elio Bava, head of the INRIM, Torino). 
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Figure 1: Phasor representation of a randomly modulated signal. 



Using the peak amplitudes instead of the rms amplitudes, we introduce the 
signal (Fig. 1 A) 



Vi{t) = Vb[l + a(t)] cos[w * + ¥>(*)] , 



(1) 



which defines the random fractional amplitude a(t) and the random phase f(t). 
The above can be rewritten as [16, 17] 



Vi (t) = Vg cos coot + x (t) cos toot — y(t) cos w^t 



(2) 



where x(i) and y(i) random processes. The two above representations of Uj(t) 
are connected by 



a{t) = V Q 



1 



Vn 



ip(t) = arg(y, Vb + x) . 



1 fractional amplitude (3) 

phase, if € (— 7r,7r] (4) 



In low-noise conditions, that is, |x/Vb| <C 1 and |y/Vb| 1, and c^(t) turn 
into 



a(i) 



a; 
To 



and <^(t) = 



y_ 
Vo 



(5) 



For calibration purposes, we need a signal with pure AM (ip — 0, thus y = 0), 
and one with pure PM (a = 0, thus x = 0), as in Fig. 1 B-C. 

Phase noise is usually described in term of S v (f), namely, the one-sided 
power spectral density (PSD) of ip(t) as a function of the Fourier frequency /. 
The physical dimension of S v (f) is rad 2 Hz _1 . Similarly, S a (f) is the PSD of 
a(t). S a (f) has the dimension of Hz -1 . The technical quantity £(/) is also 
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used, defined as £(/) = \S v (f) and given in dBc (dB carrier). The technical 
unit dBc is also used for AM noise. Trusting the stationarity and the ergodicity 
of the random process, that is, the repeatability and the reproducibility of the 
experiment, the PSD is measured as the average square modulus of the one- 
sided Fourier transform normalized for the power-type signals. The general 
background on phase noise and on frequency stability is available from numerous 
references, among which we prefer [5, 10, 27, 4]. 

1.2 State of the art 

Generally, the measurement of AM and PM noise is done with a suitable detector 
followed by a fast Fourier transform (FFT) analyzer. In most practical cases, 
the phase detector is a double balanced mixer [13, 2, 1] with the two inputs 
in quadrature. The amplitude detector [18] is a power detector, that is, a 
diode used in the quadratic region, or a double balanced mixer with the two 
inputs in phase. A balanced bridge (often referred to as 'interferometer') with 
amplification and synchronous detection of the noise sidebands is used when 
the highest sensitivity is required [24, 11]. The sensitivity is limited by the 
equivalent temperature of the instrument [9]. Improved sensitivity is obtained 
by correlation and averaging, with two separate — thus independent — systems 
that measure the same device under test (DUT) [26, 32]. The dual-bridge 
with correlation exhibits the highest reported sensitivity, limited by the thermal 
uniformity of the instrument instead of the absolute temperature [22]. The 
bridge works well for the measurement of two-port components, yet it is difficult 
to use for the oscillator PM noise, and virtually unsuitable to the oscillator AM 
noise. 

Calibration is the accurate measurement of the amplitude-to- volt age gain 
k a and of the phase-to-voltage gain k v of the detector. The calibration relics 
upon a reference of phase modulation and of amplitude modulation. A de- 
terministic modulation is often used because the uncertainty of the base-band 
instruments that follow the detector can be made small as compared to that 
of the RF/microwave section. Figure 2 shows the calibration schemes generally 
used. 

Following the scheme A, a sideband of frequency v s close to vq is added to the 
carrier in one arm of the phase detector. This is equivalent to phase modulation 
plus amplitude modulation of rms value ^P S /2P , and frequency Vb = \v s — 
P s and Pq are the sideband power and the carrier power, respectively. The 
scheme B is similar to the scheme A, but for the use of RF /microwave random 
noise instead of the sideband. The close-in noise is S a (f) = S v (f) = N/Pq, 
where N is the noise PSD. 

The schemes A and B suffer from the impossibility to divide AM from PM. 
Instead, an equal amount AM and PM is present. Hence, a standard based on 
Fig. 2A-B should be referred to as a standard of SNR (signal-to-noise-ratio), 
rather than of AM and PM noise. The problem is that the mixer used as 
the phase detector is not perfectly balanced, for the dc offset is affected by 
the input power. Hence, it turns the AM noise into near-dc noise, which is 
mistaken for PM noise. The AM noise rejection can be of 15-35 dB, depending 
on the operating power, frequency and quadrature error, on the mixer type, and 
on technology [20, 21]. Using the mixer at a sweet point of zero-sensitivity to 
AM noise, as suggested in [3, 6]. Consequently, the use of this standard in the 
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Figure 2: Generally used calibration methods. 



calibration of PM noise measurements based on the saturated mixer is incorrect. 

In the scheme of Fig. 2 C, two sinusoidal signals are sent to the mixer. These 
signals have frequency vq and v s , with a small difference vi> = \ v s — v$ \ , and the 
same power of the final application. The detected signal is a sinusoid of fre- 
quency Vb crossing the V dc axis with slope ^ — tn^k^, from which k v can be 
measured using a digital oscilloscope. Of course, this scheme can only be used 
only for the calibration of PM measurements, not for AM. The major prob- 
lem with this approach is that it is impossible to assess the mixer's unwanted 
sensitivity to AM noise. 
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The scheme of Fig. 2 D makes use of a true amplitude or phase modulator. 
This is definitely correct, provided the modulator was calibrated. Avoiding AM 
to PM crosstalk, that is, avoiding residual AM in the phase modulator and 
residual PM in the amplitude modulator, is the critical point. 

Presently, the state of the art in AM and PM noise standards is the NIST 
secondary standard of AM and PM noise [31, 14]. This standard is based on 
the scheme of Fig. 2 B, thus it makes the use of the mixer incorrect as discussed 
before. Additionally, however accurate, this standard is a secondary one. Com- 
mercial standards of this type also exist, like that of Techtrol Cyclonetics Inc. 
(TCI) [25]. 

A reference modulator (Fig. 2D) is in use at NIST 2 , calibrated with the 
secondary standard. Yet, nothing is said about how to avoid the AM to PM 
crosstalk. 

In the BIPM web site 3 , only the French LNE reports on the accuracy of phase 
noise measurements. The uncertainty is of 2 dB. However disappointing, this 
value accounts for all conditions of routine measurements. The NIST 4 reports 
about typical accuracy of 1 dB in both AM and PM noise measurements. An 
old article on the bridge [7]. The declared uncertainty is of 1.9 dB at the NPL, 
UK. 

2 Method 

This article proposes a method for the primary calibration of AM and PM 
measurements, based on the following ideas. 

• We first understand that the calibration of the detector (amplitude or 
phase) is the critical point, while the accuracy of the near-dc complemen- 
tary instruments (FFT, etc.) can be easily made negligible for all practical 
purposes. 

• A reference modulator is used to calibrate the gain of the detector. Of 
course, we need an amplitude modulator free from phase modulation, 
and a phase modulator free from amplitude modulation. The conceptual 
scheme of the modulator is shown in Fig. 3. The amplitude modulation is 
obtained by adding to the carrier Vb cos(coot) a signal x(t) cos(uJot). Thus, 
the amplitude modulation is a(t) = x(t)/Vo [Equations (1) and (2)]. Sim- 
ilarly, the phase modulation is obtained by adding a signal — y(t) cos(woi), 
so that ip(t) = y(t) /Vb 

• For our purpose, it is necessary to fix the quadrature error ip and the gain 
asymmetry e of the real I-Q modulators and detectors (Fig. 6). This can 
be done with null and differential measurements only. 

• Real phase detectors show a residual sensitivity to amplitude modulation. 
This is due to non-perfect saturation. Conversely, power detectors are 
not sensitive to phase modulation, provided the video 5 bandwidth is large 
enough to prevent the discriminator effect due to the internal memory. 
Consequently, we can use null measurement with a power detector to find 

2 Craig Nelson, Phase Noise Measurements, 2006 FCS tutorial 
3 http:/ /kedb. bipm.org/appcndixC/. . . 

4 http://ts. nist.gov/MeasuremcntServices/Calibrations/oscillators.cfm 
5 The output is usually called 'video.' 



E. Rubiola, AM-PM calibration. January 8, 2009 



RF input 

V cos(d) t) 



RF output 

V cos(to t) +xcos(o)()t) -ysin((B t) 



LO 



i-Q 

modul 



RF 



modulation sidebands 
xcos(co t) -ysin((O t) 




I-Q modulator 

LO pump 



X y modul. input 



) 




output 


1 


Q 


) • 



Figure 3: Reference modulator. 



the phase 6 of the pump. 9 is an arbitrary quantity resulting from the 
electrical layout, to be corrected in order to set the appropriate phase 
relationship between the modulation sidebands and the carrier. 

• Switching on and off the modulation sideband, the modulation is measured 
as a power ratio. Thus, ce Tms = s/P x /Po, and (p rras — WPy/Po, where Pq>, 
P x and P y are the power of the carrier, of the in-phase modulation, and 
of the quadrature modulation. 

It is to be remarked that the entire procedure proposed here relies only on null, 
differential, and ratiometric measurements. This will be discussed in detail after 
explaining some technical tools. 



3 Mixers, scalar products, and Hilbert spaces 

It is well known in the classical theory of telecommunications [33, 29, 12] that 
with the appropriate definition of the scalar product, the space of the sinusoids 
at a given frequency (a>o) has the structure of an Hilbert space. In fact 

1. the space is a inner-product space, that is, the scalar product (x,y) has 
the usual properties of conjugation and linearity, and (x,x) > 0, where 
(x, x) = only if x = 0, 

2. the norm is = <J (x, x), 

3. the space is complete. 

In such space, we can define a orthogonal unitary base, so that all signals can 
be expressed as a linear combination of the base. 

3.1 Complex signals 

The scalar product between two complex signals s(t) and r(t) is defined as 

(*,r) = i J\{t)r*{t)dt (6) 
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where T is a suitable integration time, and the superscript * denotes the complex 
conjugate. The base function is the complex sinusoid 



r(t) = e JLJot 
l|r(*)||=l 



reference 



(7) 
(8) 



Mathematically, the integration time should be T = n^, integer n, so that the 
u>o and 2o>o terms are canceled. In practice, a sufficiently long T> ^ makes 
the 2u>o terms negligible for any practical purposes. 
Let us consider the scalar product of the signal 

s(t) = Ve j " ot+v signal (9) 

on the reference r(t). Expanding Eq. (6), the scalar product v = (s,r) is 



exp(ju t + ip) exp(-jw t) 



dt 



V 



cos(u>ot + tp) + j sm(ujot + if) cos(wo^) — j sin(wo^) 



dt 



V j cos(u t) cos(<,o) — sm(Lo t) sin(^) 



+ 3 



((Jot) sin(<p) — sin(oj £) cos(^)j j 

cos(woi) — j sin(oj i)J dt 



and finally 

v = V\cos((f) + j sin(ip)~\ (after removing the 2o>o terms) 

= Vfpv 

Note that the signal s(t) can be expressed as 

s(t) = vr(t) 



(10) 
(11) 

(12) 



The important fact is that the scalar product freezes the u>o carrier frequency, 
thus it gives the phasor v — Ve^ v associated to s(t). V is the rms voltage. 

3.2 Real signals 

Dealing with real signals s(t) and r(t), the scalar product (6) turns into 

r-T 



(s,r)=^J o s(t)r(t) 



dt 



(13) 



because r(t) — r*(t), thus the complex conjugate is no longer needed. A two- 
dimension base is necessary to replace the complex sinusoid 



ri(t) — + V2cos(u)ot) 
r Q(t) = -\/2sin(w t) 



in-phase reference (14) 
quadrature reference (15) 
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The factor \[2 is necessary for ||r/|| = 1 and ||r g || = 1. Let us consider the 
scalar product of 

s(t) — vV2cos(ujQt + ip) signal (16) 

on the base functions. Expanding Eq. (13), the scalar product vi = (s,ri) is 

i r T 

vi = — 2V cos(ojot + <p) cos(ojot) dt 
1 Jo 

i f T r i 

= — / 2V cos(u) t) cos(<y9) — sm(oj t) sin(<p) cos(w i) dt 
T Jo L J 

i f T r i 

= — / 2^ cos 2 (woi) cos(^) — sin(w i) cos(wo^) cos(y) (it 

J JO L ' J 

= Vcos^) (after removing the 2w terms) 

Similarly, the scalar product vq — (s,tq) is 



V Q = ™ / 2t/cos(w t + (/?) [- sin(w i)] eft 



■T 



2V 



— sin(wot) 
dt 



cos(wo^) cos(<^) — sin(wot) sin(<£) 

1 fT f 

= — / 2V sin(w i) cos(w i) cos(y>) + sin 2 (w i) sin(<^) 
J Jo L 

= Vsin((/3) (after removing the 2uj terms) 

Note that the real signal s(t) can be expressed as 

s(t) = vin(t) +v Q r Q (t) . 



dt 



(17) 



(18) 



(19) 



Once again, the signal s(t) is completely defined by the scalar projections 
the base functions. The proof is trivial 



on 



s(t) — V cos(tp) \/2cos(a;o*) — Vsin(^) \/2sm(u t) 
= vV2cos(LU t + ip) . 



(20) 
(21) 



3.3 I-Q detector and modulator 

The I-Q detector is a hardware implementation of Eq. (17) and (18). Fig. 4 
shows the conceptual scheme of an ideal I-Q detector, based on analog multi- 
pliers. The rules of energy management and impedance matching arc still not 
considered here. In the geometrical analogy, the I-Q detector extracts the real 
and imaginary part of the input phasor. 

The I-Q modulator (Fig. 4) is similar to the detector, but the signal flow 
is reversed. In the geometrical analogy, the I-Q modulator turns vj into a real 
phasor and vq into an imaginary phasor, and builds the output phasor by adding 
the real part and imaginary part. 
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3.4 Defects of the real I-Q detectors and modulators 

Practical I-Q detectors and modulators are based on double balanced mixers 
and on baluns. See the reference [19] for a tutorial. Figures 6 and 7 show 
the defects of a practical I-Q detector and modulator, detailed underneath. In 
addition, we stress that the inputs can not be driven with a voltage source. 

DC offset and LO crosstalk. Due to the asymmetry of the diode ring or of 
the baluns, a dc offset appears at the output of the detector. The same 
phenomenon manifests as the imperfect suppression of the carrier at the 
modulator output when the inputs are set to zero. 

Loss. The mixer operates with the LO input saturated. Accordingly, the peak 
amplitude of the reference signals takes the value Vl, which generally 
differs from \J2. For a given mixer used as a detector, Vl can be calculated 
from the mixer loss £ using Vl = ^j-- The typical loss of an I-Q detector 
is I = 2.8 (9 dB, i.e., the 6 dB of a mixer plus the 3 dB of the input power 
splitter), we calculate Vl = 700 mV. 



Gain asymmetry. The mixer loss is not the same for the I and Q arm. A 
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Figure 7: Errors of the real I-Q detector. 

difference of 5% is common. 

Orthogonality error. The phase shift of the pump in the Q arm deviates by 
an angle tf> from the quadrature. This error, which can be of a few degrees, 
is frequency dependent. 

Frame orientation error. The arbitrary length of the cable connecting the 
LO pump to the reference turns into an arbitrary phase 9, which affects 
both the I and Q arms. 



Correcting these defect is a central point in our work. 
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Figure 9: Measurement of the orthogonality and symmetry defect of a real I-Q 
detector. 

4 Correction of the I-Q detector 

This section describes how to compensate the defects of a I-Q detector, turning 
it into a virtually ideal one, as shown in Fig. 8. The offset is removed by adding 
an appropriate dc. The the orthogonality defect ip and the gain asymmetry e 
are fixed by taking a linear combination of the two outputs. The phase of the 
LO pump is still let arbitrary, as it results from the electrical layout. 



4.1 Measurement of the errors i\y (orthogonality) and e 
(symmetry) 

First, the dc offset is compensated for after inspection with a dc voltmeter, 
with the RF input terminated. Afterwards, the matrix D corrects for ip and e. 
The method to obtain the elements of D from the measurements is shown in 
Figures 9, and described underneath. 
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We feed a low-power sideband 

v r (t) = cos(w + ^b)t (22) 
= cos(woi) cos(u>bt) — sin(wot) sin^t) (23) 

in the RF input of the detector. For the sake of simplicity, the amplitude 
coefficients are omitted. The frequency uo s = ut, + ojo is chosen for w& to be an 
appropriate audio frequency (1-2 kHz). The LO input is saturated by a signal 
cos(w i) of suitable power. The output signals 

wj(t) = *cos(w 6 i) (24) 
VQ(t) = ^(1 + e) [sm(ujbt) cos(ip) — cos(u)bt) sin(^>)] (25) 

= i(l + e)sin(w 6 i- V) (26) 

are measured with a lock-in amplifier, which gives tp and e. The lock-in takes 
the / signal as the oscillator reference. Most lock-in amplifiers measure only the 
amplitude of the input, not the amplitude of the oscillator reference. If so, input 
and oscillator reference are to be interchanged in order to get both amplitudes. 
Interchanging the lock-in inputs also helps to get an accurate measure of ip. 



4.2 Correction of the orthogonality and symmetry errors 

The detector is corrected by introducing the matrix D at the output (Fig. 8). In 
the remaining part of this section we denote with v\ (t) and v'q (t) the outputs of 
the I-Q detector, and with vi(t) and vq(z) the outputs of the corrected detector, 
after the matrix D. Thus, 



VI 




dn 


di2 




'v'j 






d 2 i 


d 2 2_ 







(27) 



The elements of D are calculated by separating the cos(wo^) and sin(wot) parts 
of the input signal, i.e., the real and imaginary parts of the input phasor. Thus, 
when the RF input is 



v r (t) — cos(wq^) cos(w(,i) (real part) 



(28) 



the output 



= 2 cos M 

v' Q (t) 



— - ( 1 + e) cos (ujbt ) sin ip 



(29) 
(30) 



must be transformed into 



v i(t) = ^cos(uj b t) 
v Q (t) 



2 




(31) 
(32) 
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This is equivalent to state that 



d 2 i d 2 2 



\ cos(u) b t) 
— \ ( 1 + e) cos (u>bt) sin 





cos(u>bt) 








holds for all t, which yields 



d n = 2 
d 12 = 
d2i - d 22 (l + e) simp = . 



(1 + e) sin(w^) cos ip 



Similarly, when the RF input is 

v r {t) = — sin(o>oi) sin(wbt) 

the output 

v'i(t) = 
v' Q (t) 

must be transformed into 

w/(t) = 

v Q(t) = 2 s'M^bt) ■ 
This is equivalent to state that 



\ (1 + e) s\n(tjjbt) cos tp 



(imaginary part) 



dn 


d\2 




d 2 i 


d22_ 












sin(wbt) 



holds for all t, which yields 



d 2 2 = 



(1 + e) cost/ 1 
Finally, combining the above results we get 



D 



1 

sin-0 



(33) 



(34) 
(35) 
(36) 



(37) 

(38) 
(39) 

(40) 
(41) 

(42) 
(43) 



(44) 



_costp (1 + e) costp. 
4.3 Correction of the pump phase 

As a result of the previous sections, the I-Q detector is corrected for the ampli- 
tude asymmetry e and for the quadrature error ip. The phase 9 of the LO pump 
is still arbitrary. This can be corrected (Fig. 10) by introducing at the matrix 



R 



cos f 
sin 8 



— sin b 
cos 8 



(45) 
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Figure 10: Correction of the arbitrary phase 9 of the LO pump. 



which rotates the output by 9. This requires an absolute reference of angle, 
which we still have not. For the sake of completeness, we show the mathematical 
derivation of the matrix R. 

We start by feeding the signal 



v r (t) = xcos(uj()t) — ysm(ujQt) 



(46) 



at the RF input. The I-Q internal pump signals are cos(uJot+9) and cos(wn,i+#)- 
As before, we denote with v'j(t) and Ug(t) the outputs of the I-Q detector (now 
corrected by the matrix D), and with Vi(t) and i>g(f) the outputs of the corrected 
detector, after the matrix R. In this section we discard the factor 2 introduced 
by the matrix D, as we had an ideal I-Q detector without correction. Thus, the 
output signals are 



v'j(t) = [xcos(u>ot) — ysm(u>ot)\ cos(ujot + 9) 

= \x cos(a>o^) — y sin(wot)] [cos(wo£) cos 9 — sin(woi) sin 9] 

= — [xcos(9) + ysin6>] (after removing the 2ojo terms) 

Vg(t) — [cccos(woi) — V sin(wo^)] {— sin(ujQt + #)} 

= [xcos(woi) — y sin(woi)] {— [sin(woi) cos + cos(woi) sin (9] } 



— - [a;sin(r5) — ycost 



(after removing the 2wq terms) 



(47) 
(48) 

(49) 

(50) 
(51) 

(52) 



This is a rotation of the vector x, y by —9, which is corrected with the transfor- 
mation 



(53) 
(54) 



VI 




rn 


ri2 




v/ 












- V Q. 






cos t 


) - 


sin 8 






sin 6 


cost? 





Wi] 







Hence we have found Eq. (45), QED. 
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Figure 11: Correction of the real I-Q modulator. 



4.4 Accuracy considerations 

The above described process moves the measurement uncertainty from the 
RF/microwave section to the near-dc section, where the instrument uncertainty 
can be made negligible. Additionally, we remark the following facts. 

• Modern lock-in amplifiers measure the real and imaginary part of the input 
signal after in-phase and quadrature down-conversion to dc. Hence, the 
measurement of the small quadrature defect ip is actually a null measure- 
ment. This is definitely true if the process is reiterated for best accuracy. 

• Using a digital lock-in, the lock-in internal quadrature reference relies on 
sampling, and ultimately on the internal clock. This is free from the phase 
error of an analog detector. 

• The measurement of the small gain error e is a differential measurement, 
which is limited by noise, rather than by the calibration of the lock-in. 

5 Correction of the I-Q modulator 

This section describes how to compensate the defects of a I-Q modulator, turning 
it into a virtually ideal one. The phase of the LO pump is still let arbitrary, as 
it results from the electrical layout. The method, shown in Fig. 11, is similar to 
that used to fix the I-Q detector. 

5.1 Parameter measurement 

First, the LO crosstalk is eliminated by adding an appropriate dc offset, after 
inspection with a spectrum analyzer. Afterwards, the matrix M corrects for the 
orthogonality defect ip and for the gain asymmetry e. The method to obtain 
the elements of M from the measurements is shown in Figure 12, and described 
underneath. 

Feeding a signal 



x = V m cos(ui m t) 

y = o 



(55) 
(56) 
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Figure 12: Measurementof the orthogonality and symmetry defect of a real I- 
modulator. 



at the modulator inputs, we calculate 9 after measuring the output signals 



v x = -V m cos(u; m t) cos I 



v y = -V m coa(cj m t) sin( 



Similarly, feeding a signal 



x = 

V = V m cos(uj m t) 



(57) 
(58) 



(59) 
(60) 



at the modulator inputs, we measure the output signals 



V m (l + e ) cos(w m t) sin(# + ijj) 



'J y = ^V m (l + e) cos(w m t) cos(6» + ip) 



(61) 
(62) 



from which we calculate e and ip. 

It is to be remarked that the use of small signals of the form V m cos(w m i) 
allows the accurate measurement of v x and v y with a lock-in amplifier. 

A variant of this method consists of inserting a phase shifter in the RF line 
between the modulator and the detector, and to adjust the phase for 8 = 0. 
After eliminating 9 with a null measurement, Eq. (61)-(62) provide a more 
accurate estimation of e and ip. The phase shifter should be introduced in 
the RF line, rather than in the LO pump of the modulator or of the detector, 
because these devices have a residual sensitivity to the pump power, and the 
phase shifter introduces a loss. 
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5.2 Fixing the defects 



The system, as shown in Fig. 12 is equivalent to a channel described by the 
relationship 



v x 




Cll 


Cl2 




X 


Vy_ 




C21 


C22_ 




y. 



where the matrix C is 



C = 



icos6» -§sin(6> + V) 
\sin9 \{l + e)cos(6> + tp) 



(63) 



(64) 



Hence, the system is corrected by inserting the matrix C 1 (the inverse of C) 
at the input 



c- 1 



2cos(9 + tp) 2 sin(6» + ?/>)' 
cos tp cos tp 

2sin# 2cosfl 

(l + e)cosi/> (1 + e) cos tp 



(65) 



cos ip 



cos(6 + -0) sin(6» + tp) 
sin cos 9 



l + e 



l + e 



(66) 



so that C _1 C = / (identity matrix). As we have still not defined the measure- 
ment layout, there is no point in correcting 8. Thus we set 9 = in C . The 
matrix C _1 |e = o is denoted with M 



M 



cos tp 
1 



(67) 



(1 + e) cos tp 

In summary, the matrix M fixes e and tp, and lets the pump phase 6 arbitrary. 



6 Building the reference AM-PM modulator 

As a result of Section 5 we have a virtually ideal I-Q modulator to be inserted in 
the reference AM-PM modulator (Fig. 3). Yet, the electrical layout introduces 
constant arbitrary phase shifts, which we represent as a single phase 9 at the 
LO input of the I-Q modulator. In this section we measure and correct 9. 



6.1 Parameter measurement 

We first observe that phase modulators have a residual sensitivity to amplitude 
modulation and noise, while the power detector is insensitive to phase modu- 
lation and noise, provided the video bandwidth is large. Therefore, we use a 
power detector to set 9 for the I-Q modulator to add pure phase modulation to 
the carrier when it receives a signal at the Q input. This id detected as a null of 
amplitude modulation at the detector output. Having corrected the geometry 
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Figure 13: Setting the phase of the modulator pump. 



errors of the I-Q modulator, the channel / adds pure amplitude modulation to 
the carrier. 

Fig. 13 shows the method to measure 9. The I-Q pump is cos(wo^ + S), for 
the I-Q RF output signal is vj cos(ojo^ + 6) — V Q sin(woi + &). This is added to 
the carrier Vocos(u>ot). Setting the switch to /, we get 



vi(t) = V m cos(uj m t) and v Q (t) = 
at the input of the I-Q modulator, and 

Vi(t) = V COs(wrjt) + V m COs(lV m t) Vq COs(w t + 9) 

at the detector input. The detected signal is 

Mt) = \ V Q + \ v rl + VoV m cos(6) cos(w m i) + cos{2u m t) 
Similarly, setting the switch to Q we get 

vi(t) = and y(i) = vq cos(u) m t) 
at the input of the I-Q modulator, and 

Vi(t) = Vq cos(uJot) — V m cos(uj m t) Vq sin(w £ + 9) 
at the detector input. The detected signal is 

Mt) = \Vq + \vl - V Vm sin(0) cos(w m t) + hft cos(2u; m t) 



(68) 
(69) 

(70) 

(71) 
(72) 

(73) 



The lock-in amplifier selects the u> m term in Equations (70) and (73), from which 
we calculate 9. 

6.2 Fixing the defects 

The result of the measurement process is summarized as follows 



VI 




'Vm 










v d = cV V m cos( 



(74) 
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Figure 14: Setting the phase of the modulator pump. 



AM PM 



and 







' " 






v m 



Vd 



-cVoVm sin 9 



which can be rewritten as the scalar product 

j 



Vd 



Vi 
Vq 



[cos 9 — sin 0] 



vi 
V Q. 



(75) 

(76) 
(77) 



where the amplitude coefficients cVoV m are dropped in order to simplify the 
formalism. 

The arbitrary phase shift 9 is compensated for by introducing a matrix R at 
the input of the I-Q modulator (Fig. 14) 



VI 




ru 


ri2 






VQ_ 




r n 


ri2_ 




Vy_ 



so that 



Vd = [c 



rn 


ri2 




V x 




ri2_ 




Vy_ 



(78) 



(79) 



The matrix R has to satisfy two conditions: when the input is v x = 1, v y = 
the output is a modulation in phase with the carrier, which is detected 



cost 



- sin 



8) 



ru 


ri2 




1" 


r n 


r i2_ 








1 



(AM) 



(80) 



and when the input is v x = 0, v y — 1 the output is a modulation in quadrature 
with the carrier, which is not seen by the power detector 



[cos 9 — sin #] 



ru 


ri2 




"0" 


fu 


ri2_ 




1 



= 



(PM) 



(81) 
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The 



solution is 



R = 



cos 



e 



sm 



e 

6 



(82) 



— sm 



(9 



cos 



7 Calibration checklist 

The following list summarizes the sequence of actions for proper calibration. 

1. Correct the I-Q detector (e and ip) 4. 

2. Correct the I-Q modulator (e and ip) 5. 

3. Build the reference AM-PM modulator 6. 

• chose the appropriate levels. 

• fix 9. 

• measure the modulation depth. 

4. Optional, needed for AM: fix 8 of the detector. 

5. Use the reference modulator to calibrate a measurement, system 

8 Modulation depth and accuracy 

First we address a question left aside, why calibrating the system using a sinu- 
soidal modulation instead of random or pseudorandom noise. The reasons are 
detailed unerneath. 

• Quadrature modulation is not true PM, it produces some AM. Yet, the 
residual AM occurs at the frequency 2iv m , which is rejected by a lock-in. 

• The noise bandwidth can be made narrow. 

• Can reject the background noise of the system. 

Though all these features can also be achieved with pseudorandom modulation, 
the sinusoidal modulation gives them in a simple and straightforward way. 

Figure 15 shows a variant of the reference modulator (first introduced in 
Fig. 3), adapted to the measurement of the modulation depth. The modulator 
is used in the following ways. 

Pass through. This is used to measure the noise of a DUT. For highest accu- 
racy of the signal path, the modulator is not removed. This makes sense 
only if the modulator loss in series to the DUT can be tolerated. 

Calibrate. Normal operation of the reference modulator, which provides a cal- 
ibrated amplitude or phase modulation. 

Self calibration. The reference modulator is calibrated in two steps, as de- 
tailed in Section 8.2. 
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Figure 15: Measurement of the modulation depth. 



8.1 Phase and frequency modulation 

Sinusoidal phase and frequency modulations can be used as a reference of angle 
modulation. Unfortunately, in (virtually) all practical cases this method does 
not provide simultaneously the reference unmodulated signal together with the 
modulated signal. 

The basic phase modulation is 

v(t) = cos [w t + tp(t)] (basic PM) . (83) 

In the presence of a sinusoidal modulation ip(t) = msm(uj m t) of modulation 
index m, the signal writes 

v(t) = cos [(Jot + msm(ui m t)] (sinusoidal PM) , (84) 

which can be expanded as 

v(t) = cos(ujot) cos [msm(a> m i)] — sm(uj Q t) sin [msin(w m t)] (85) 

= Jo(m) cos(wo^) ~ 2J\{m) sin(w m t) sin(uj t) + . . . (86) 

= Jo(m) cos(wo^) + Ji(m) sin(o'o + Um)t — J\{m) sin(o;o — w m )t (87) 

~ cos(wo^) + ^msin(wo + oj m )t — |msin(wo — w m )t (small m) . (88) 

Expanding the sinusoidal PM, we have used the following properties 

cos [z sin 6] — Jo(z) cos 8 + . . . and sin [zsin#] = 2Ji(z) sinO + . . . , 

sin(a) sin(6) = | [cos(a — b) — cos(a + b)] , 

Jo(z) = 1 + . . . and Ji(z) = \z + . . . for m — > . 

Finally, equating the phase modulation tp(t) = msin(u) m t) to the phase tp(t) = 
J(Aw)(t) dt obtained from the frequency modulation yields 

(Aw)(t) = (Acj)j, cos(w m t) (sinusoidal FM) , (89) 
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and therefore 



m = = (modulation index) . (90) 

8.1.1 Commercial synthesizers 

Using commercial synthesizers, the frequency modulation is preferred to the 
phase modulation because of easier calibration. The DDS is a remarkable ex- 
ception to this statement because the phase modulation can be obtained nu- 
merically by adding a term to the accumulator content, before addressing the 
look-up table. 

The modulation index can be measured precisely measuring (Av) p in static 
conditions with frequency counter, by switching a dc input and trusting the 
linearity of the VCO for small frequency deviation. 

Another method, known as the Bessel null method, consists of increasing m 
until the carrier nulls. From (86), the carrier nulls at the first zero of Jo (to), 
which occurs at m = 2.405. Observing v(t) with a spectrum analyzer, the carrier 
drops suddenly in a sharp interval. 

A third method consists of the direct measurement of m from the ratio of the 
carrier power P divided by the sideband power P s , using a spectrum analyzer. 
At small modulation index, (88) gives 

8.2 Measurement of the modulation depth 

We measure the carrier power and the modulation power separately. A reference 
attenuation i is introduced in the carrier path, so that the carrier power Po/i 2 
and the DSB power are of the same order of magnitude. This ensures precise 
power ratio measurements even with commercial power meters. The modulation 
is 

Fp~ rp^ 

frms = \ or "rms = \ ~ET ( 92 ) 

V -H) V -f"o 

where P x applies to AM, P y to PM. The modulation power (P x or P y ) ac- 
counts for both sidebands. A detailed explanation of the method, restrained 
to the phase modulation, is given underneath. The extension to the amplitude 
modulation is trivial. 

The quadrature-modulated signal can be written as 

v(t) = V cos(w t) - V s sin[(w + u b )t] - V s sin[(w - u b )t] (93) 
= Vb cos(uj n t) + 2V S sin(w i) cos(wf,t) (94) 

The phase modulation (p(t) associated to the above signal is 



ip(t) = arctan 



2^ cos(w b t) 



(95) 



2^ cos(ui b t) linearized for small V s /V (96) 
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thus 

^rms = yF& (97) 
V 2 

The carrier power is Pn = wrh, where i? is the load resistance. Similarly, the 
power of each sideband is P s = thus P y = 2;^ for both sidebands. Hence, 
the power ratio ^ is 



Py 2Ro _ VJ 



' -2^x^-2^- (98) 



By virtue of Eq. (97), the above is equal to (</? 2 ), which proves Equation (92). 



8.3 Mixer response and linearity 

The sole mixer of our interest in this article is the double-balanced mixer based 
on a ring of Schottky diodes (see [19] for a tutorial). The mixer used as an 
analog modulator is expected to respond linearly 

v RF = ai IF (99) 

where a is an appropriate coefficient. Actual mixers deviate from linearity, and 
saturate at an output power of some 1 mW (0 dBm) . We measured a mixer Mini 
Circuits ZFM-2 driven at the IF input by a current generator in the following 
conditions: 

LO frequency 100 MHz 
LO power 5 mW (7 dBm) 

RF load 50 fl 

IF resistance 1 kf2 

Due to the exponential behavior of the diodes and to the even symmetry, 
similar to that of a differential pair, it is natural to model the mixer as 

vrf — a\ tanh(a2i/F) pure tanh(x) model (100) 



or as 



vrf — o-i tanh(o 2 i/F) + a ^iiF tanh(.x) model with dissipation (101) 

The second model differs from the first by the insertion of a series resistance, 
which accounts for the diode and transformer resistance. The experimental data 
and the least-square approximations are shown in Figure 16. For the sake of 
completeness, we also try a polynomial fit (Fig. 17). Table 2 summarizes the 
coefficients of the least-square fits. 

Looking at the above models, we notice that the tanh(x) model is superior to 
the polynomial fit, and that adding the dissipative term has a small beneficial 
effect. Table 1 gives an idea of the error Avrf/vrf if the the full Equation 
(100) or (101) is replaced with the tangent. 

We expect that the model (100) or (101) is a factor 10 more accurate than 
the linear approximation, that is, 10~ 3 at ijp = 1 mA and 10 -5 at iiF = 100 
yuA. We expect uncertainty and errors of the same order when the mixer is used 
as a synchronous detector. Of course, this is the accuracy of the model, which 
does not include the uncertainty of the instruments used to measure the mixer. 
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Figure 16: Double-balanced mixer used as an analog modulator. Top: experi- 
mental data, polynomial fit, and small-signal (linear) approximation. Bottom: 
discrepancy between model and linear approximation. 



Table 1: Linearity defect [Eq. (101)] of a double-balanced mixer used as a 
modulator. 



ilF &V RF /v R f 


vrf Prf 


0.1 10" 4 
0.316 10~ 3 
1 10~ 2 


3.91 0.305 -35.2 
12.4 3.05 -25.2 
39.1 30.5 -15.2 


mA (dimcnsionlcss) 


mVtms dBm 



8.4 Power levels and signal-to-noise ratio (SNR) 

The choice of power levels depends on the application, and chiefly on the power 
of the signal to be calibrated. Figure 18 (derived from Fig. 15) provides a first 
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Figure 17: Double-balanced mixer used as an analog modulator. Top: ex- 
perimental data, 3rd and 5th degree polynomial fit, and small-signal (linear) 
approximation. Bottom: discrepancy between model and linear approximation. 



example, in which the DUT is suitable to a power of 10 mW (+10 dBm) and 
shows an attenuation of 1 dB. The directional coupler in which the modulation is 
injected shows an attenuation of 1 dB (intrinsic and dissipative). Thus Pq = 6.3 
mW (+8 dBm). We measure Pq by inserting a reference 40 dB attenuator (round 
value), so that P /l 2 = 800 nW (-31 dBm). In this We try i IF = 316 /iA as 
a compromise between linearity and SNR. In this conditions the modulation 
power is P y = 316 nW, after 10 dB attenuation of the directional coupler. 
The equivalent noise power at the amplifier input is Pm — FkTB. Assuming 
a noise figure F — 1.25 (1 dB) and a noise bandwidth B = 1 MHz, we get 
Pn = 5xl0 -15 W, thus a signal-to-noise ratio 

SNR = = 3 - 05xl0 ' 7 r = 6.05x10 s (87.8 dB) (102) 
Pat 5.04x10- 15 v ; y ' 
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Table 2: Coefficients of the numerical fits for a ZFM-2 mixer used as a modulator 
(see Fig. 16 and 17). 

polynomial approx, RMS voltage 
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-1.444e5 
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a7 
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al a3 a5 
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Figure 18: Tentative estimation of signal levels and SNR. 



It is to be remarked that B impacts only on the SNR, thus on the assessment 
of the uncertainty, not on the modulation depth. A narrower B can be used, 
provided the bandpass region is flat enough to attenuate the carrier and the 
sidebands equally. Table 3 compares some calibration scenarios. 

8.5 Tentative estimation of the uncertainty 

Though the method relics only on null and ratio measurements, having primary- 
metrology facility for power and impedance on site helps in achieving the highest 
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Table 3: Expected SNR of the scheme shown in Fig. 18. 



%l F 


Prf 


SNR 


0.1 


0.305 -35.2 


7.79xl0 3 


77.8 


0.316 


3.05 -25.2 


2.46xl0 4 


87.8 


1 


30.5 -15.2 


7.79xl0 4 


97.8 


mA 


/iW dBm 


(dimcnsionlcss) 


dB 



Table 4: Tentative estimation of the uncertainty, based on commercially avail- 
able component and instruments. The values given below arc thought to be 
conservative. 



parameter and conditions 


value 




power ratio measurement 


11.6X10" 3 


(0.1 dB) 


(commercial power meter) 






RF path 

(couplers, cables etc.) 


23xl0" 3 


(0.2 dB) 


reference attenuator 


5.8xl0- 3 


(0.05 dB) 


(Weinschel, 40 dB) 




mixer and detector linearity 


l.OxlO" 3 




Section 8.3 






null measurements 


l.OxlO" 3 




(commercial lock-in, 10 bit) 






signal-to-noise ratio 


l.OxlO" 3 




Section 8.4 






worst case total 


43.6xl0" 3 


(0.37 dB) 


rms total 


26.5xl0- 3 


(0.23 dB) 



accuracy. We attempt to estimate the uncertainty under the hypothesis that all 
the calibration is based on commercial instruments and that no primary metrol- 
ogy is available on site, which is our case. Table 4 summarizes the uncertainties, 
keeping on conservative values. 

9 Application to the bridge (interferometric) in- 
struments 

At the end of the conceptual work I realized that my bridge (interferometric) 
instruments for the measurement of AM and PM noise [23] contain almost all 
the blocks needed for the calibration. Therefore, the calibration process can be 
added with a minimum of changes. This concept is shown in Figure 19, where 
the new blocks (in color) are added to the instrument. 
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Figure 19: The calibration method fits well into the bridge (interferometric) 
noise measurement instruments. The notation is that of the reference [23], 
where the dual bridge is published. In color, the changes needed for calibration. 



10 Where these ideas come from 

Though at superficial sight our scheme has some similarity to a scheme patented 
by F. L. Walls at NIST [30], theory of operation, calibration procedure, and use 
are radically different. The key idea of synthesizing a modulated signal by 
adding carrier and sidebands is so old that I can't track back the original idea. 
I wish to conclude this Section with a warning against patents. The SI units, 
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by definition, can not be patented. In the official documents of La convention 
du metre and its child BIPM it is written clearly that the SI units are defined 
and designed for them to be free and accessible. Some readers could remember 
that one of the major objection against the definition of the second based on 
the Earth revolution around the Sun, in vigor in 1964-68, was that it was not 
sufficiently accessible. 
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